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ABSTRACT 



In this paper I present three new results of astronomical interest concerning the 
theory of Abel inversion. 1) I show that in the case of a spatial emissivity that is 
constant on toroidal surfaces and projected along the symmetry axis perpendicu- 
lar to the torus' equatorial plane, it is possible to invert the projection integral. 
From the surface (i.e. projected) brightness profile one then formally recovers the 
original spatial distribution as a function of the toroidal radius. 2) By applying the 
above-described inversion formula, I show that if the projected profile is described 
by a truncated off-center gaussian, the functional form of the related spatial emis- 
sivity is very simple and - most important - nowhere negative for any value of 
the gaussian parameters, a property which is not guaranteed - in general - by 
Abel inversion. 3) Finally, I show how a generic multimodal centrally symmet- 
ric brightness distribution can be deprojected using a sum of truncated off-center 
gaussians, recovering the spatial emissivity as a sum of nowhere negative toroidal 
distributions. 

keywords: Methods: analytical - Methods: numerical - Methods: data analysis 
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1. INTRODUCTION 



A common problem in astronomy is the deprojection of a given surface bright- 
ness distribution. Unfortunately this problem is degenerate, i.e., different spatial 
emissivities can originate the same surface (projected) brightness distribution. As 
a consequence, any inversion procedure is invariably based on a (more or less) 
arbitrary choice of the underlying geometry of the spatial emissivity. A help to 
this unpleasant situation comes from two guide rules: in choosing the geometry of 
the spatial emissivity one uses symmetry properties (if any) of the brightness dis- 
tribution, and after deprojection discards the assumed geometry if this produces 
a somewhere negative spatial emissivity. 

One of the simplest cases is given by a surface brightness characterized by 
central symmetry, i.e., described by a function I(R) of the projected radius R. The 
natural assumption is a spherically symmetric spatial emissivity, £ = £{r) > 0, 
where r is the spherical radius. It is a well known result that in this case the 
projection operator and its deprojection are given by an Abel integral equation: 

I{R) = 2 r g(r)rdr 
1 ' J R Vr 2 -R 2 K ' 

and 

'n dl(R) dR _ f(r t ) 
J r dR - r 2 y/r?-r 2 \ ' 

where r t is the truncation radius, i.e., I(R) = for r > rt, and for untruncated 
distributions r t — > oo. It is important to note that if I(r t ) > the recovered 
emissivity is (weakly) divergent at the edge of the sphere, due to the second therm 
on the r.h.s. of eq. (2). A fundamental problem of Abel inversion is that the 
positivity of the recovered distribution is not guaranteed. As an example of this 
fact, one can consider I(R) to be a gaussian or an off-center gaussian: while in the 
first case the deprojected emissivity - being still a gaussian - is everywhere posi- 
tive (see, for applications, Bendinelli, Ciotti, & Parmeggiani 1993, and references 
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£(r) = — 
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therein), in the second case it can be easily shown that £{r) diverges negatively 
for r — > 0. 

What other symmetries for £ are compatible with a surface brightness I = 
I(R) and are still of astrophysical interest? The natural generalization of the 
spherical symmetry is the cylindrical one. Many astronomical objects are char- 
acterized by this geometry: lenticular and spiral galaxies, accretion disks, plasma 
tori, planetary rings, some planetary nebula etc.. In this paper we are interested 
in particular in the problem of deprojecting a I(R) which presents some off-center 
maxima, and so the natural assumption about the emissivity distribution £ is the 
toroidal one. 

In Section 2 I show that in the case of a toroidal emissivity distribution pro- 
jected along its symmetry axis it is possible to generalize the inversion formula 
(2). Then I show that in the particular case of a projected brightness described 
by a truncated off-center gaussian, £ has an extremely simple functional form, 
and at variance with the spherically symmetric case, it is nowhere negative for 
any choice of the gaussian parameters. In Section 3 I propose the use of the found 
physically admissible IS pair to obtain a (finite) series expansion of the emis- 
sivity for a generic centrally symmetric surface brightness profile. In a following 
paper (Bendinelli et al. 1999), an application of this method to the Planetary 
Nebula (PN) A 13 is described. This PN is just one case among many others of 
astrophysical interest to which the method produces useful results on the object's 
structure. 

2. THE INVERSION FORMULA FOR TOROIDAL SYMMETRY 

Let us start by obtaining the analogous formula of eq. (1) for the projection along 
the symmetry axis of an emissivity distribution stratified on toroidal surfaces. The 
natural coordinates are the cylindrical, (R, cp, z) , where z = is the torus equa- 
torial plane. The relation between cylindrical and toroidal coordinates (r, <p, #) is 
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given by 



R = R + r sin i?, <p = (p, z = r cos 



(3) 



(see Figure 1). 

Figure 1 T7ie describing parameters of toroidal geometry, and their relations used 
for the projection and deprojection. 

We assume independence* of £ from # and cp, and so each isoemissivity surface is 



placed in the equatorial plane. Moreover, for the sake of generality let us assume 
that the emissivity is non-zero only for < r < rt < Ro- For rt = Ro we have a 
so-called full torus. With the previous assumptions, I(R) > for |_R — i? | < rt 
and I(R) = outside. From very simple geometric arguments one obtains the 
analogous formula of eq. (1): 



With respect to the variable R, the l.h.s. of eq. (4) is a generalized Abel integral, 
but it is not directly invertible, since (R — R ) 2 is not strictly increasing with R 
(see Gorenfio & Vessella 1991, p. 24). Since the function's profile is symmetric with 
respect to R D , we can use a branch of I (i.e., we use R > R D ) and we define a new 
variable s := R — R > with I+(s) := I(R + s). In this way, the integral is 
formally identical to the projection operator in eq. (1) and can be inverted: 



* The request of independence of 8 from <p is unnecessary, and the following 
discussion can be easily generalized to distributions I = I(R,cp), obtaining £ = 



labeled by its toroidal radius r = ^/ z 2 + (R — R ) 2 around the circle of radius R Q 







£{r,cp). 
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Thus, we proved that assuming the surface brightness distribution to be the pro- 
jection on its equatorial plane of a toroidally stratified emissivity, it is formally 
possible to invert the projection integral, recovering the emissivity as a function of 
the toroidal radius. 

3. TRUNCATED OFF-CENTER GAUSSIANS 

As an application of eq. (5) we invert a surface brightness profile described by an 
off-center gaussian truncated at r t : 



I(R) = S \ exp 



(R-R ) 2 ( r 



2(7 2 



2 



exp(-^)|>, \R-R \<r t . (6) 



Note that i+(r t ) = 0, and so the unpleasant divergence at the edge of the torus 
is avoided. The total luminosity L associated to / is given by its surface integral 
over the annulus R Q — rt < R < R Q + r t : 



L = SR D a (2tt) 3 / 2 



Erf fe)"^ exp ("^ 



(7) 



where Erf (a:) = 'i/^/lx exp(— t 2 )dt is the error function. The formula for the 
luminosity density given using eqs. (5)-(6) results to be: 



W = ^-exp -— Erf \ . (8) 




The deprojection formula can be verified by inserting eq. (8) in eq. (4) and 
then evaluating the integral. Note that the spatial emissivity given by eq. (8) is 
everywhere positive, for any choice of the gaussian parameters (S, R , r t , a). One 
can then conclude with the following analogy: off-centered gaussians in toroidal 
symmetry correspond to centered gaussians in spherical symmetry. 

Having found a nowhere negative IS pair that satisfies the Abel inversion for 
a particular toroidal density distributions, the successive step is the natural exten- 
sion of a previous work (Bendinelli et al. 1993), where a multigaussian expansion 
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of observed profiles is described, under the assumption of spherical symmetry: here 
it is sufficient to remember that by using gaussian functions one avoids the com- 
putational difficulty of direct Abel numerical inversion, which belongs to the class 
of unstable inverse problems (Gorenflo & Vessella f99f ; Craig & Brown 1986). 

So one can assume that a given profile I(R) is expanded as a sum of truncated 
off-center gaussians with different parameters (S, R , rt, <r)i for i = 1,...,N. The 
parameters are easily computable using the Newton-Gauss regularized method 
which is a powerful iterative non linear fitting technique (Bendinelli et al. 1987; 
Bendinelli 1991). Clearly the parameters are constrained to satisfy the request 
that the integral over all the projection plane of the fitted brightness distribution 
must be equal to the original one, i.e. the total luminosity of the system L = 
J^Lj, where Lj are given by eq. (7). The associated spatial emissivity is then 
approximated by the sum of N distributions as in eq. (8) . The developed technique 
will be applied to the deprojection of the galactic PN A 13, which is characterized 
by a well defined ring-shaped morphology (Bendinelli et al. 1999). 

4. CONCLUSIONS 

The results of this work are the following: 

1) It is shown that it is possible to extend the classical Abel inversion for spheri- 
cally symmetric systems to toroidal density distributions projected along the torus 
symmetry axis. 

2) It is also shown that an off-center truncated gaussian function gives rise to a 
well behaved spatial emissivity, i.e., the emissivity is a very simple function of the 
toroidal radius. More important, the emissivity is non-negative over all the space 
for any choice of its parameters, a property not guaranteed by the Abel inversion. 

3) Finally, it is proposed the use of the found IS pair for the recovering of the 
spatial emissivity of any centrally symmetric projected distribution as a sum of 
toroidal density distribution, after a non-linear parameter fitting. 
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